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ABSTRACT 

One of the puzzles of the dark energy problem is the (first) cosmological coincidence problem, 
namely, why does our universe begin the accelerated expansion recently? why are we living in an 
epoch in which the dark energy density and the dust matter energy density are comparable? On 
the other hand, cosmological observations hint that the equation-of-state parameter (EoS) of dark 
energy crossed the phantom divide Wde = — 1 in the near past. Many dark energy models whose 
EoS can cross the phantom divide have been proposed. However, to our knowledge, these models 
with crossing the phantom divide only provide the possibility that Wde can cross —1. They do not 
answer another question, namely, why crossing the phantom divide occurs recently? Since in many 
existing models whose EoS can cross the phantom divide, Wde undulates around —1 randomly, why 
are we living in an epoch Wde < —1? This can be regarded as the second cosmological coincidence 
problem. In this work, the cosmological evolution of the vector-like dark energy interacting with 
background perfect fluid is investigated. We find that the first and second cosmological coincidence 
problems can be alleviated at the same time in this scenario. 
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I. INTRODUCTION 

Dark energy problem flj has been one of the most active fields in modern cosmology, since the discovery 
of accelerated expansion of our universe ^^i^i^^ |61j ■ In the observational cosmology, the equation- 
of-state parameter (EoS) of dark energy Wde = Pde/ Pde plays a central role, where pde and pde are its 
pressure and energy density, respectively. To accelerate the expansion, the EoS of dark energy must satisfy 
Wde < —1/3. The simplest candidate of the dark energy is a tiny positive time-independent cosmological 
constant A, whose EoS is —1. However, it is difficult to understand why the cosmological constant is 
about 120 orders of magnitude smaller than its natural expectation, i.e. the Planck energy density. This is 
the so-called cosmological constant problem. Another puzzle of the dark energy is the (first) cosmological 
coincidence problem, namely, why does our universe begin the accelerated expansion recently? why are we 
living in an epoch in which the dark energy density and the dust matter energy density are comparable ? This 
problem becomes very serious especially for the cosmological constant as the dark energy candidate. The 
cosmological constant remains unchanged while the energy densities of dust matter and radiation decrease 
rapidly with the expansion of our universe. Thus, it is necessary to make some fine-tunings. In order 
to give a reasonable interpretation to the (first) cosmological coincidence problem, many dynamical dark 
energy models have been proposed as alternatives to the cosmological constant, such as quintessence 0,0, 
phantom 9, 10, jj, k-essence ,12.] etc. 

Recently, by fitting the SNe la data, marginal evidence for Wde{z) < — 1 at redshift z < 0.2 has been 
found ^3 ■ In addition, many best-fits of the present value of Wde are less than — 1 in various data fittings 
with different parameterizations (see 0| for a recent review). The present observational data seem to 
slightly favor an evolving dark energy with Wde crossing —1 from above to below in the near past |l5j |. 
Obviously, the EoS of dark energy Wde cannot cross the so-called phantom divide Wde — for quintessence 
or phantom alone. Although it seems possible for some variants of k-essence to give a promising solution 
to cross the phantom divide, a no-go theorem, shown in shatters this kind of hopes. In fact, it is not 
a trivial task to build dark e nergy m odel whose EoS can cross the phantom divide. To this end, a lot of 

efforts 111 El El El 111 IMM El IS El 111 111 113 ^ have been made. However, to our 

knowledge, many of these models only provide the possibility that Wde can cross —1. They do not answer 
another question, namely, why crossing phantom divide occurs recently? Since in many existing models 
whose EoS can cross the phantom divide, Wde undulates around —1 randomly, why are we living in an 
epoch Wde < — 1 ? It can be regarded as the second cosmological coincidence problern [331 . 

The so-called second cosmological coincidence problem was seriously discussed in |23|_for the first time. 
In 33], the key point is the trigger mechanism, similar to the case of hybrid inflation |34| and the model 
by Gong and Kim [ssl ] . In the hybrid dark energy model considered in js^ , a quintessence and a phantom 
are employed. The feature of spontaneous symmetry breaking plays a critical role. In the first stage, the 
phantom field a is trapped to cr = 0, while the quintessence field slowly rolls down for a long time. In 
this stage, the effective EoS of hybrid dark energy remains larger than —1. When (/) reaches a critical value 
(pc, a phase transition takes place and the phantom field a is triggered to climb up its potential. And then, 
the quintessence field is trapped to (fi — 0. Crossing the phantom divide occurs between the moment 
of phase transition (j) — (pc and the moment of </> being trapped at </> = eventually. The effective EoS 
remains smaller than —1 in the major part of the second stage. When a continuously climbs up, it will 
eventually arrive at a critical value Oc- The quintessence field 4> is triggered to roll down again. In the 
third stage, whether the EoS remains smaller than —1 or changes to be larger than —1 depends on the 
profile of the effective potential of hybrid dark energy in this stage. Thus, the avoidance of the big rip is 
possible for suitable model parameters. We refer to the original paper for more details. (Also, see |60ll 
for other discussion on the second cosmological coincidence problem, soon after our previous paper [23 ■) 

Although the scalar field is used extensively because of its simplicity, the vector field also gets its 
applications in modern cosmology. Here, we only mention a few, such as in [s^ l37| . the vector field is 
considered as the source which drives the inflation; in [ssL l39l ] the Lorentz-violated vector field and its 
effects to universe are studied; the quantum fluctuations of vector fields produced at the first stage of 
reheating after inflation is also studied in ^3]; cosmology of massive vector fields with SO (3) global 
symmetry is investigated in ^1]; and also see m 23] and references therein for the literature on the 
magnetic fields in the universe. It is found that non-linear electromagnetic field can drive the acceleration 
of the universe 0] . The vector field can be a good dark matter candidate reproducing flat rotation curves 
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in dark halos of spiral galaxies (see j43 for example) . Of course, the vector field is also a viable dark energy 
candidate |4y, |43 ■ The effects of vector field dark energy candidate on the cosmic microwave background 
radiation and the large scale structure are discussed in |58j |. 

In this work, the cosmological evolution of the vector-like dark energy proposed in |4(il | interacting with 
background perfect fluid is investigated. In fact, in the cosmological evolution of the vector-like dark 
energy with inverse power-law potential and without interaction with background matter was studied, by 
means of directly numerical and approximate solutions of the equation of motion in two limits of matter 
domination and vector-like dark energy domination. Differing from "4^, we investigate the cosmological 
evolution of the vector-like dark energy and background perfect fluid in this work by means of dynamical 
system |48|. And we study the models with not only inverse power- law potential but also exponential 
potential. Furthermore, we consider the cases of vector-like dark energy interacting with background 
perfect fluid, while the interaction terms are taken to be four different forms which are familiar in the 
literature. We find that the first and second cosmological coincidence problems can be alleviated at the 
same time in this scenario. 

This paper is organized as follows. In Sec. ^ we will briefly present the main points of the vector- 
like dark energy model proposed in [4^ . In Sec. IIIII we give out the equations of the dynamical system 
of vector-like dark energy with interaction to background perfect fluid for the most general case. That 
is, we leave the potential of vector-like dark energy and the interaction form undetermined. We will 
investigate the dynamical system for the models with inverse power-law and exponential potentials in 
Sec. IIVI and respectively. In each case with different potential, we consider four different interaction 
forms between vector-like dark energy and background perfect fluid. The interaction forms are taken to be 
the most familiar interaction ones considered in the literature. In Sec. lVIl the first and second cosmological 
coincidence problems are discussed. Finally, brief conclusion and discussion are given in Sec. IVIII 

We use the units h ^ c = 1 and = SttG throughout this paper. 



II. 



VECTOR-LIKE DARK ENERGY 



In fact, the "vector-like dark energy" is a so-called "cosmic triad" (see |43), which is a set of three 
identical vectors pointing in mutually orthogonal directions, in order to avoid violations of isotropy. Fol- 
lowing |46j | , we consider the case of vector-like dark energy minimally coupled to gravity, and the action is 
given by 




3 r 



a=l 



(1) 



where g is the determinant of the metric g^^, R is the Ricci scalar, Fj^^ = d^A^ — = g'^^A'^A^, 

and Cm is the Lagrangian density of matter fields ip. Latin indices label the different fields (a, 6, . . . = 1,2,3) 
and Greek indices label different spacetime components (/i, j^, . . . = 0,1,2,3). Actually, the number of 
vector fields is dictated by the number of spatial dimensions and the requirement of isotropy j4^. The 
Latin indices are raised and lowered with the flat "metric" 6ab- The potential term V (A"^) explicitly 
violates gauge invariance. In 46], it is argued that this kind of "cosmic triad" can naturally arise from 
gauge theory with a single SU{2) gauge group. 

From action one can get the energy-momentum tensor of the cosmic triad and the equations of 
motion for the vectors A"; as 



(A) J. 

J- II. 



E 



pa pap 

l-Lp V 



9-^5^ /la Aa 



-^F^p^F^P" + V {A 



a2\ 



(2) 



and 



dV 



(3) 
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respectively. We consider a spatially flat Friedmann-Robertson- Walker (FRW) universe with metric 

ds^ ^ ~dt^ + a^{t)dx.'^, (4) 

where a is the scale factor. An ansatz for the vectors, which can be compatible with homogeneity and 
isotropy, is 

A\ = S\A{t) ■ a. (5) 

Thus, the three vectors point in mutually orthogonal spatial directions, and share the same time-dependent 
length, i.e. A"'^ = A'^A'"'' = A'^{t). Substituting Eqs. © and JH into Eq. Q, one obtains 

A + ZHA+[h' + ^)a+'^^0. (6) 

where H = d/a is the Hubble parameter, and a dot denotes the derivative with respect to the cosmic time 
t. The Friedmann equation and Raychaudhuri equation are given by, respectively, 

^2 ^2 

H'^ = -jPtot = y (pa + Pm) , (7) 

and 

^2 ^2 
H ^ — {ptot + Plot) = ^ (pa + Prn +PA+ Pm) , (8) 

where Pm and pm are the pressure and energy density of background matter, respectively. The energy 
density and (isotropic) pressure of the vector-like dark energy are given by 

PA^l[A + HAy +W{A^), (9) 

1 / . \ 2 dV 

Pa = ^{a + HA) -3V{A')+2^A^ (10) 

respectively. Noting that d/a = H^ + H, one can check that Eq. ^ is equivalent to the energy conservation 
equation of vector-like dark energy, namely pA + 'iH{pA + Pa) — 0. 

The most remarkable feature of the vector-like dark energy is that its EoS wa = Pa/ Pa can be smaller 
than —1, while possessing a conventional positive kinetic term. This is thanks to the additional term in 
proportion to dV/dA^ in pA- While the energy density pA is positive, one can find that the condition for 
WA < ~1 is 

Thus, dV/dA^ < is necessary . There are other interesting issues concerning the vector-like dark 

energy. We refer to the original paper for more details. 



III. DYNAMICAL SYSTEM OF VECTOR-LIKE DARK ENERGY INTERACTING WITH 

BACKGROUND PERFECT FLUID 

As mentioned in Sec. ^ in this work we will generalize the original vector-like dark energy model [4^ 
to include the interaction between the vector- like dark energy and background matter. The background 
matter is described by a perfect fluid with barotropic equation of state 

Pm = WmPm = (7 " l)Pm, (12) 
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where the barotropic index 7 is a constant and satisfies < 7 < 2. In particular, 7=1 and 4/3 correspond 
to dust matter and radiation, respectively. We assume the vector-like dark energy and background matter 
interact through an interaction term C, according to 



PA + (pa + pa) = -C, 

Pm + 3H {prn + Pm) = C, 



(13) 
(14) 



which preserves the total energy conservation equation ptot + 3i? {ptot + Ptot) = 0. It is worth noting that 
the equation of motion ^ should be changed when C ^ 0, a new term due to C will appear in its right 
hand side. 

Following 1^ I51I IH^ , we introduce following dimensionless variables 



kA 



kVv 
7^' 



'^V Pm 



kA 



(15) 



By the help of Eqs. (|7|)- (|10|l . the evolution equations (|13|l and H14|l can then be rewritten as a dynamical 
system |4^, i.e. 



x' ^6 
y' 62/ 

z' = 6z 



{x + uf + ^z^ + e {x + u)~ 2eu-^ ~3x-2u-Ci, 



{x + uY + — + ( 1 H — xy '^u 



{x + uY 
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1 



e 



e- 



C2 



where 



and 



e 



dV 
IPdA^' 



K^C 



(x + u) ^, 



C2 



zC 



18H^ - 2Hp,,_ 

a prime denotes derivative with respect to the so-called e-folding time A'^ = Ina, and we have used 



H 



= 6 



{x + uf + ^z^ + Q 
4 



The Friedmann constraint Eq. Q becomes 



(x + u) + jY 



+ z' 



The fractional energy densities of vector-like dark energy and background matter are given by 



3if2 



[x + uY + lY 



fir 



2,m 



respectively. The EoS of vector-like dark energy and the effective EoS of the whole system are 



WA 



PA 

PA 



{x + uY - 3y2 + 46 



(x + u) + 



and 



^e// = = ^AWa + ^rnWm, 

Ptot 



(16) 

(17) 

(18) 
(19) 

(20) 
(21) 

(22) 
(23) 
(24) 

(25) 
(26) 
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respectively. From Eq. it is easy to see that the condition for wa < — 1 is 



{x + uf + Q<Q, (27) 

which is equivalent to Eq. Hll() . Finally, it is worth noting that y > and z >Qhy definition, and in what 
follows, we only consider the case of expanding universe with H > 0. 

It is easy to see that Eqs. ((TB|) - (tT^ become an autonomous system when the potential V (A^) is chosen 
to be an inverse power-law or exponential potential and the interaction term C is chosen to be a suitable 
form. Indeed, we will consider the model with an inverse power-law and exponential potential in Sec. IIVI 
and0 respectively. In each model with different potential, we consider four cases with different interaction 
forms between vector-like dark energy and background perfect fluid, which are taken as the most familiar 
interaction terms extensively considered in the literature H3j IEj> 1^ lEE lE^ i namely 

Case (I) C = 3aHpm, 

Case (II) C = SpHptot = SpH {pA + Pm) , 

Case (III) C = r]Kp„iA, 

Case (IV) C = 'iaHpA, 

where a, /3, rj and a are dimensionless constants. 

In the next two sections, we firstly obtain the critical points (x, y, z, u) of the autonomous system by 
imposing the conditions x' ~ y' = z' = u' = 0. Of course, they are subject to the Friedmann constraint, 



namely 3 



{x + u) +y 



+ z = 1. We then discuss the existence and stability of these critical points. An 



attractor is one of the stable critical points of the autonomous system. 



IV. MODEL WITH INVERSE POWER-LAW POTENTIAL 



In this section, we consider the vector-like dark energy model with an inverse power-law potential 



V{A^) =Vo{k'AY\ (28) 
where n is a positive dimensionless constant (required by the condition dV/dA^ < 0). In this case, 

e = -^y^. (29) 

One can obtain the critical points {x,y, z,u) of the dynamical system Eqs. (|16|I - H19|I with Eqs. (|28|l and 
H29(l by imposing the conditions x' — y' — z' = u' = 0. Note that these critical points must satisfy the 
Friedmann constraint H23|l . y > 0, z > and the requirement of x, y, z, u all being real. To study the 
stability of these critical points, we substitute linear perturbations x ^ x + Sx, y y + Sy, z ^ z + 5z, 
and u ^ u + 6u about the critical point {x,y,z,u) into dynamical system Eqs. (|16|I - H19|I with Eqs. (|28|l 
and 1)29(1 and linearize them. Because of the Friedmann constraint 123(1 . there are only three independent 
evolution equations, i.e. 



Sx' 



9{n + 2){x + uf - 2nu~^{x + u) + [ + n] z 



Sx 



+2z 



7 -1- n I {x + u) 



n , 
-3(1 



-2\ --2 
Z \ u 



2nxu 



Sz- 



3n- 2 



(18 + 9n + nu~'^) {x + uf + Q7 



Su — SCi, 



Sz' = &{n + 2)z{x + u)Sx + 



3(n + 2){x + uf + [-l + n] {iz^ - l) 



+6(n + 2)z{x + u)Su + SC2, 
Su' = Sx, 



(30) 



(31) 
(32) 



where SCi and SC2 are the linear perturbations coming from Ci and C2, respectively. The three eigenvalues 
of the coefficient matrix of the above equations determine the stability of the corresponding critical point. 
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A. Case (I) C = SaHpm 

In this case, Ci = ^ {x+u)~^ and C2 = |o;z. The physically reasonable critical points of the dynamical 
system Eqs. ()16|) -()19 | with Eqs. (|28|l and (|29|l are summarized in Table|l| Next, we consider the stability of 
these critical points. Substituting (5Ci = —^z^{x + u)~^Sx + az{x + u)~^Sz — ^z^{x + u)~'^6u, 5C2 — |a<5z, 
and the corresponding critical point {x,y,z,u) into Eqs. H3() |l - H32|l . we find that Point (P.I.I) is always 
unstable; Point (P.I. 2) exists and is stable under condition a < 7; Point (P.I. 3) is unstable if it can exist. 
The unique late time attractor (P.I. 2) has 

= 1, = 0, WA = -1, Weff = -1, (33) 

which is a vector-like dark energy dominated solution. 



Label 


Critical Point (a;, y, z, u) 


P.I.I 


0, 0, 0, ±^ 


P.I.2 


0' ]/ 3{2+n)' 0' 3(2 + n) 


P.I.3 


0, 0, ^1-3^^, ±^37^ 



TABLE I: Critical points for Case (I) C — SaHpm in the model with inverse power-law potential. 



Case (II) C = SpHptot = 3I3H {pA + p„ 



In this case, Ci = ^{x + u) ^ and C2 — ^Pz ^. We present the physically reasonable critical points of 
the dynamical system Eqs. (fT?))l - ltT??|l with Eqs. I2HJ and (|^ in Table ITU where 



12/3 



37 



(34) 



-4(i + w)-2 (Sx + Su), SC2 



-(3z '^Sz, and the corresponding critical point 



Then, we substitute SCi 

{x,y,z,u) into Eqs. (|30|I - H32|) to study its stability. We find that Point (P.II.l) exists and is stable under 

- 1 1 and 7 > 8/3, which is out of the range < 7 < 2; Point (P.II.2) 



conditions < (3 < ^— ^ 



12 



4-1-37 



3(37-4) 



exists and is stable under conditions (3 < min |o, 

and is stable in a proper parameter-space [H^ . 
The late time attractor (P.II.2) has 



4-37 
12 



(4+37)" 



9(4-37)^ 



1 



I and 7 < 4/3; Point (P.II.3) exists 



an = ^(1+ri), 



iVeff^-[il + n){3j-i) + 2] 



(35) 



which is a scaling solution. The late time attractor (P.II.3) has 
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WA = -I- 



Pi 

7-/3' 



= -1, 



(36) 



which is a scaling solution also. Note that Weff of attractor (P.II.2) is larger than —1, since < 7 < 2 
and < ri < 1 which is required by Eq. (|34|l and its corresponding fim, while wa of attractor (P.II.3) is 
smaller than —1, since < /? < 7 is required by its corresponding 17™. 
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Label 



Critical Point {x, y, z, u) 



P.II.l 
P.II.2 
P.II.3 



0, 0, [i(l-.,)]''',±[|(l + -i)]''' 

47+)9(-4+37) fil 



6(2+n)7 



2„(7-/3)-3/37 
6{2+n)7 



TABLE II: Critical points for Case (II) C = SPHptot 
potential, ri is given in Eq. 13411 . 



3l3H {pA + Pm) in the model with inverse power-law 



C. Case (III) C = -qt^PmA 



In this case, Ci = "^(2; + u) ^xz'^ and C2 = ^^rjxz. The physically reasonable critical points of the 
dynamical system Eqs. l(TO|l - lfT^ with Eqs. if^ and if^ are shown in Table ^] Substituting 5Ci = 



-^(x + u) [1 - (a: + u) ^x] (5a;-|- 4/|?7(a; + u) ^xzJ^; 



rjx5z, 



and the corresponding critical point (x,?/, z, u) into Eqs. H3()|l --I|^^2 |l . we find that Point (P.III.l) is always 
unstable, while Point (P.III.2) is always stable. The unique late time attractor (P.III.2) has 



VIa = 1, ^rn = 0, WA 

which is a vector-like dark energy dominated solution. 



-1, 



-1, 



(37) 



Label 



P.III.l 
P.III.2 



Critical Point (a::, y, z, u) 



0, 0, 0, ±^ 



3{2+n) ' 



0, ± 



3(2+n) 



TABLE III: Critical points for Case (III) C = rjKpmA in the model with inverse power-law potential. 



D. 



Case (IV) C = ZaHpA 



In this case, Ci = fcr [(a; + uY + y^] [x + uy^ = ^{l- z'^){x + uy^ and C2 = fcr [{x + u 
|(T (z^^ — zY The physically reasonable critical points of the dynamical system Eqs. H16(I - (I19|I with 
Eqs. if^ and are shown in Table |^ where 



r2 



-37 



(38) 



Next, we consider the stability of these critical points. Substituting 5Ci 



az{x + u)-^6z - |(x -h uY"^ (1 - z^) Su, SC2 



-^{x + u)-^ {l-z^)6x 



+ 1) 6z, and the corresponding critical point 



{x,y,z,u) into Eqs. H3U|) - (|32|I . we find that Point (P.IV.l) is unstable if it can exist, while Point (P.IV.2) 
exists and is stable in a proper parameter-space |57l |. The unique late time attractor (P.IV.2) has 



^A 



7 + CT 



WA 



-1, 



(39) 



which is a scaling solution. Note that wa of attractor (P.IV.2) is smaller than —1, since cr > is required 
by its corresponding flm- 
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Label 



P.IV.l 
P.IV.2 



Critical Point {x, y, z, u) 



0, 0, (3r2)l/^ ±(i -ra) 

± 



1/2 



7(4+3tT) 



6{2+n)(7 + o-) ' y 7 + CT ' 



7(2n-3<7) 



6(2 + n){7 + CT) 



TABLE IV: Critical points for Case (IV) C = 3aHpA in the model with inverse power-law potential. r2 is given 
in Eq. (EHJ. 



V. MODEL WITH EXPONENTIAL POTENTIAL 

In this section, we consider the vector-like dark energy model with an exponential potential 

V (A^) = Vo exp {-Xk^A^) , (40) 
where A is a positive dimensionless constant (required by the condition dV/dA'^ < 0). In this case, 



e 



-3Xu^y^ 



(41) 



One can obtain the critical points {x,y,z,u) of the dynamical system Eqs. (|1()|I - H19|I with Eqs. I|4()|l and 
H41(l by imposing the conditions x' — y' ~ z' ^ u' = 0. Note that these critical points must satisfy the 
Friedmann constraint H23|l . y > 0, z > and the requirement of x, y, z, u all being real. To study the 
stability of these critical points, we substitute linear perturbations x ^ x + 5x, y y + Sy, z ^ z + 6z, 
and u ^ u + Su about the critical point {x,y,z,u) into dynamical system Eqs. (|16|I - H19|I with Eqs. H4()(l 
and 1)41(1 and linearize them. Because of the Friedmann constraint H23|) . there are only three independent 
evolution equations, namely 



Sx' = 6 



(1 + 3Xu^) {x + u)2 - 2Au(x + u)+ (Au^ + ^^z^- Xu^ - ^ 
3{x + u) (Aw^ + - Au Sz + 2| 18Xu{x + u)^ + 3 {9Xu^ + 3 - A) (x + u)^ 



dz' 



-6Xu{x + u)+ 6Xu{x + u)+3 (Au^ + |) - A (z^ - l) + -7 - l| Su - 5Ci 

12z (1 -I- 3Xu^) {x + u)Sx -t- 6 (1-1- 3Xu^) {x + uf + [Xii^ + |) (32^ - l) 
+\2z [3Xu{x + uf + (1 + 3Am^) (x + u)+ Xu (z^ - l)] 5u + SC2, 



Sz 



Su' ~ 5x, 



(42) 



(43) 
(44) 



where SCi and 6C2 are the linear perturbations coming from Ci and C2, respectively. The three eigenvalues 
of the coefficient matrix of the above equations determine the stability of the corresponding critical point. 



A. Case (I) C = SaHpm 

In this case, Ci = %{x+u)^^ z"^ and C2 = Jaz. The physically reasonable critical points of the dynamical 
system Eqs. ((16|l - (|19|l with Eqs. ((4()|) and l(4l)l are summarized in TableFVI Next, we consider the stability of 
these critical points. Substituting (5Ci = —^z-^{x + u)~^Sx + az{x + u)~^Sz — ^z'^{x + u)~'^Su, SC2 = ^aSz, 
and the corresponding critical point {x,y,z,u) into Eqs. (|42|l -144 () . we find that Point (E.I.I) is always 
unstable; Point (E.I. 2) exists and is stable under conditions a < 7 and A > 1; Point (E.I. 3) is unstable if 
it can exist. The unique late time attractor (E.I. 2) has 

= 1, r2,„ ==0, WA^ -1, Weff = -1, (45) 

which is a vector-like dark energy dominated solution. 
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Label 


Critical Point {x, y, z, u) 


E.I.I 
E.I.2 
E.I.3 


0, 0, 0, ±^ 
0, 0, ^l-^,±^^ 



TABLE V: Critical points for Case (I) C = 3aHpm in the model with exponential potential. 



B. 



Case (II) C = SpHptot = 3/3// (pa + Pm) 



In this case, Ci — ^{x + u) ^ and C2 = ^. We present the physically reasonable critical points of 
the dynamical system Eqs. (fT?)|l - l(T^ with Eqs. and (|IT|l in Table IVII where 



= V-3l3j^X + (7 + /3A - jXy 



(46) 



Then, we substitute 5Ci ^ —^{x + u) ^ {6x + Su), SC2 = —^P^ '^8z, and the corresponding critical point 
{x,y,z,u) into Eqs. H42|) - (|44|l to study its stability. We find that Point (E.II.l) exists and is stable under 



conditions < f3 < 



4-37 
12 



4+37 



3(37-4) 



- 1 > and 7 > 8/3, which is out of the range < 7 < 2; Point (E.II.2) 



in{0, ^ 



37 



exists and is stable under conditions /3 < niin 

(E.II.4) exist and are stable in a proper parameter-space 
The late time attractor (E.II.2) has 



(4+37)" 



1 



9(4-37)^ 

jSTj . respectively. 



I and 7 < 4/3; Points (E.II.3) and 



^e// = g[(l + ri)(37-4) + 2] 



(47) 



which is a scaling solution. The late time attractors (E.II.3) and (E.II.4) both have 



WA 



1 - 



/37 
7-/3' 



"eff 



(48) 



which are scaling solutions also. Note that Weff of attractor (E.II.2) is larger than —1, since < 7 < 2 
and < ri < 1 which is required by Eq. H34|) and its corresponding il™, while wa of attractors (E.II.3) 
and (E.II.4) are both smaller than —1, since < /3 < 7 is required by their corresponding Qrn- 



Label 



Critical Point {x, y, z, u) 



E.II.l 
E.II.2 
E.II.3 
E.II.4 



0, 0, [l(l-rO]^/^±[|(l + n)]^/' 
0, 0, [i(l + n)]^/^±[i(l-rO]^/' 



0- V^z^^z^' 

^' 67A ' 



7 ' 



7(-l + A)-/3A + r3 
67A 

7(-l + A)-/3A-r3 
67A 



TABLE VI: Critical points for Case (II) C = 3f3Hptot ~ 3/3// {pA + pm) in the model with exponential potential, 
ri and rs are given in Eqs. 1341 and 1461 . respectively. 
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C. Case (III) C = rjupmA 



In this case, Ci = -^{x + u) ^xz^ and Ci ~ ^^rjxz. The physicaUy reasonable critical points of 
the dynamical system Eqs. ((T?)|l - l(T?^ with Eqs. (|^ and are shown in Table IVIII Substituting 



5Ci = -^{x + uY'^z'^ [l - (i + u)~^x\ 5x + \l'l-n(x + u)''^xz5z - -^(x + u)-^xz'^Su, 6C; 



rjz5x + 



rjxSz, and the corresponding critical point {x,y,z,u) into Eqs. (|42() - (|44|l . we find that Point (E.III.l) 

is always unstable, while Point (E.III.2) exists and is stable under condition A > 1. The unique late time 
attractor (E.III.2) has 



flA = 1, ^ni= 0, WA 

which is a vector-like dark energy dominated solution. 



Weff 



(49) 



Label 


Critical Point {x, y, z, u) 


E.III.l 
E.III.2 


0, 0, 0, ±^ 



TABLE VII: Critical points for Case (III) C = rjKpmA in the model with exponential potential. 



D. Case (IV) C = SaHpA 

In this case, Ci = |ct [{x + + y^] [x + uy^ = |(1 - z2)(a; + u)"i and C2 = fcr [{x + uf + y^] ^-i = 
|(T — zY The physically reasonable critical points of the dynamical system Eqs. (|16f) -()19 |l with 
Eqs. I@n|) and (|^ are shown in Table IVIlH where 

Ti = v/-37Act(7 + ct) + (7 - 7A + (t)2. (50) 

Then, we consider the stability of these critical points. Substituting 5Ci — — ^(x+u)"^ (l — z^) 5x—az{x+ 
u)~^5z ~ ^ {x + u)^^ (1 — z^) Su, 6C2 — ^^<^ (^^^ + 1) 6z, and the corresponding critical point {x, y, z, u) 
into Eqs. we find that Point (E.IV.l) is unstable if it can exist, while Points (E.IV.2) and (E.IV.3) 

exist and are stable in a proper parameter-space (57.] . respectively. The late time attractors (E.IV.2) and 
(E.IV.3) both have 

^ (J 

flA = ; , = ■ , WA = -l-(J, Weff^-l, (51) 

7 -I- CT 7 + CT 

which are scaling solutions. Note that wa of attractors (E.IV.2) and (E.IV.3) are both smaller than — 1, 
since ct > is required by their corresponding flm- 



VI. THE FIRST AND SECOND COSMOLOGICAL COINCIDENCE PROBLEMS 

As is well known, the most frequently used approach to alleviate the (first) cosmological coincidence 
problem is the scaling attractor (s) in the interacting dark energy scenario (see [iM Isollsills^. Isslls^. Isslls^ 
for examples). The dark energy can exchange energy with the background matter (usually the cold dark 
matter), through the interaction between them. The most desirable feature of dynamical system is that 
the whole system will eventually evolve to its attractors, having nothing to do with the initial conditions. 
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Label 


Critical Point {x, y, z, u) 


E.IV.l 
E.IV.2 
E.IV.3 


0, 0, (3r2)l/^±(i-r2)^/^ 

p /7(l + A) + <T-r4 1 /7(-l + A)-<T + r4 

V 6A(7+ct) ' Y 7 + <t' V 6A(7+<t) 
Q /7(l + A) + <T + r4 / 1 /7(-l + A)-a-r4 
^' Y 6A(7+ct) ' Y 7 + -^Y + 



TABLE VIII: Critical points for Case (IV) C = 3aHpA in the model with exponential potential. r2 and r4 are 
given in Eqs. H38|l and J^OJi respectively. 

Therefore, fine-tunings are needless. When the universe is attracted into the scaling attractor, a balance 
can be achieved, thanks to the interaction. In the scaling attractor, the effective densities of dark energy 
and background matter decrease in the same manner with the expansion of our universe, and the ratio 
of dark energy and background matter becomes a constant. So, it is not strange that we are living in an 
epoch when the densities of dark energy and matter are comparable. In this sense, the (first) cosmological 
coincidence problem is alleviated (see |48l Wl Isil IH^ IH^ 154 .55. .56] for examples). 

On the other hand, if the scaling attractor also has the property that its EoS of dark energy is smaller 
than —1, the second cosmological coincidence problem is alleviated at the same time. However, this is 
impossible in the interacting quintessence or k-essence scenario. Although the attractor's EoS is smaller 
than —1 in the interacting scalar phantom scenario, it is impossible to cross the phantom divide Wde — —I, 
since the EoS of scalar phantom is always smaller than —1. Fortunately, the EoS of vector- like dark energy 
can be smaller than —1, while possessing a conventional positive kinetic term |4^ . in contrast to the 
scalar phantom. Of course, the EoS of vector- like dark energy can be larg er than —1 also. Thus, crossing 
the phantom divide is possible in the vector-like dark energy model |4y|. As is explicitly shown in this 
work, for suitable interaction forms [for instance. Cases (II) and (IV)], regardless of the model with inverse 
power-law or exponential potential, there are some attractors with wa < — 1 while their corresponding 
Qa and f2m are comparable in the interacting vector- like dark energy model. In the Case (IV), all stable 
attractors have these desirable properties. Even in the Case (II), we can choose the model parameters 
to avoid the attractor with wa > —1, and the scaling attractor(s) with wa < —1 becomes the unique 
late time attractor(s). So, for a fairly wide range of initial conditions with wa > —1, the universe will 
eventually evolve to the scaling attractor (s) with wa < — 1. Similar to the ordinary way to alleviate the 
(first) cosmological coincidence problem, the second cosmological coincidence problem is alleviated at the 
same time. 



VII. CONCLUSION AND DISCUSSION 

In summary, the cosmological evolution of the vector-like dark energy interacting with background 
perfect fluid is investigated in this work. We find that the first and second cosmological coincidence 
problems can be alleviated at the same time in this scenario. Our results obtained here may support the 
vector field to be one of the viable dark energy candidates. In particular, the feature of the vector-like 
dark energy that its EoS wa = Pa/ PA can be smaller than —1 while possessing a conventional positive 
kinetic term is very attractive. While considering the interaction between the vector-like dark energy and 
background matter, the first and second cosmological coincidence problems can be alleviated at the same 
time. This is a profitable support to the vector-like dark energy. Of course, there are many remaining 
works to make this scenario more concrete, especially to fit the observational data to determine the realistic 
model parameters, which is beyond the main aim of the present work. 

The other issue is concerning the fate of our universe. It is easy to see that for all cases considered in this 
work, all stable attractors have Weff > — 1. Although the EoS of vector-like dark energy can be smaller 
than —1, the big rip never appears in this model. This is also in contrast to the ordinary phantom- like 
models. 
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Finally, we would like to mention that the gauge invariance is violated for the potential forms V (A"^) 
taken in this paper. Although in |46j| . it is argued that this kind of "cosmic triad" can naturally arise from 
gauge theory with a single SU (2) gauge group because of the equations of motion can be written down 
in a gauge invariant form, one should be careful to this problem. For the cases of inverse power-law and 
exponential potential, however, it does no longer hold, which means a particular gauge has been taken. 
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APPENDIX 

The particular parameter-space for the existence and stability of these critical points is considerably 
involved and verbose. Since our main aim here is just to point out the fact that it can exist and is stable, 
we do not present those very long expressions for the corresponding parameter-space. Of course, one can 
easily work out with the help of Mathematica. 

Instead of giving concrete expressions, we here give just some examples to support our statement. For 
the case of inverse power-law potential fSec. lIVIl . Point (P.II.3) of Case (II) exists and is stable for parame- 
ters 7 = 1, /3 = 1/3, and n = 2, with eigenvalues {-0.866239, -1.36688-^3.29093, -1.36688 + ^3.29093}, 
while it is unstable for parameters 7 = 1, /? = 1/3, and n — 1, with eigenvalues {—1.26235, 2.13117 — 
i 2.04255, 2.13117 i 2.04255}; Point (P.IV.2) of Case (IV) exists and is stable for parameters 7 = 1, 
cr = 1/2, and n = 2, with eigenvalues {-1.70481, -1.6976 + i 2.60705, -1.6976 - i 2.60705}, and for 
parameters 7 = 1, cr = 1/2, and n = 3, with eigenvalues {-1.1769, -2.32821 -h i 2.93245, -2.32821 - 
i 2.93245}, while it is unstable for parameters 7 = 1, cr = 1/2, and n = 1, with eigenvalues 
{2.15106 - i 1.12309, -2.80212, 2.15106 + z 1.12309}. Similarly, for the case of exponential potential 
(Sec. Point (E.II.3) of Case (II) exists and is stable for parameters 7=1, /3 = 1/3, and A = 6, 
while Point (E.II.4) is unstable for the same parameters; Point (E.IV.2) of Case (IV) exists and is stable 
for parameters 7 = 1, cr = 1/2, and A = 5, while Point (E.IV.3) is unstable for the same parameters. Note 
that in the above examples, we chose the demonstrative parameter 7 = 1 for dust matter, and /3 = 1/3 or 
cr = 1/2 to make flm = 1/3, which are more realistic. 
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